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theory corresponding to minimal walking technicolor. We find evidence that the pion de- 
cay constant vanishes in the infinite volume limit and that the dependence of the chiral 
condensate on quark mass m q is inconsistent with spontaneous symmetry breaking. These 
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I. MINIMAL CONFORMAL GAUGE THEORIES 



Depending on the number of flavors, matter representation and colors, non-abelian 
gauge theories are expected to exist in a number of distinct phases, classifiable accord- 
ing to the force felt between two static sources. The knowledge of this phase diagram 
is relevant for the construction of extensions of the Standard Model (SM) that in- 
voke dynamical electroweak symmetry breaking llS] • It is & l so useful in providing 
ultraviolet completions of unparticle models {1,14, 5|. 

"Minimal walking technicolor" and similar models employ fermions in higher di- 
mensional representations of the new gauge group 

la LUl- It is thought that 

some of these theories will develop a non-trivial infrared fixed point (IRFP) for a 
small number of flavors [s, The presence of a bona fide IRFP requires the van- 
ishing of the beta function for a certain value of the coupling. However, it may be 
possible (at least in perturbation theory) to find a scheme where the beta function 
has a zero yet no IRFP actually exists; indeed there are known examples in super- 
symmetric theories, when the beta function is written in 't Hooft's scheme 121 ] . On 
the other hand, if the beta function is written in a scheme that uniquely and cor- 
rectly determines scheme-independent quantities at the fixed point — such as the 
anomalous dimension (scaling exponent) of the fermion mass operator — then it is 
a "physical" beta function. We discuss such a beta function in this article — the 
conjectured all-orders beta function (cf. Eq. (12.21) below). It vanishes at g = such 
that Po — ^T(r)Nfj(g^) = where 7(g 2 ) is the anomalous dimension of ipip and T{r) 
is the Dynkin index of the representation r and Nf is the number of flavors. Since 
the one-loop coeffient (3$ is universal, it can be seen that in this scheme the vanishing 
of the beta function leads to an unambiguous result for j(gl), which is physical. 

Historically a nearly conformal behavior has been identified with the slow rise of 
the coupling constant, in an unspecified renormahzation scheme, astheenergy scale 
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It can be 



is reduced. Such a slow rise was termed walking behavior |l3l . 

mm 

shown that this is a scheme dependent statement (a nice illustration was made in 
17||). We expect, however, that a large anomalous dimension in an on-shell scheme is 
meaningful in the desired fashion — it generates a condensate that is large compared 
to the scale u Aq C d" of the theory. 

In order to better establish the location of the conformal window in minimal 



walking technico 
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or models there have recently been a number of lattice stud- 
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24. 
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26| . These investigations indicate that these gauge 



Searches for the conformal window in theories with fundamental representation quarks have also 
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theories are nearly or actually conformal, as predicted in 0, [ll|. Conformality at 
large distances implies scale invariance, which forbids a chiral condensate. Stable 
lattice simulations require a small but nonzero fermion mass, which ensures that the 
lattice theory will possess a condensate. The condensate vanishes when the mass 
is extrapolated to zero at finite volume, so lattice measurements of the condensate 
must be analyzed carefully to disentangle the infinite volume, zero mass continuum 
behavior from the effects of small but non-zero quark masses and finite volumes. In 
principle, we should be able to distinguish true scale invariance from walking behav- 
ior via a careful study of the spectrum and chiral condensate. The point is that the 
two cases will show different behavior as the mass and inverse volume are sent to 
zero, where the latter extrapolation should be performed first. This will be the aim 
of the current study which focuses on the minimal walking technicolor model. This 
is an 577(2) gauge theory with two Dirac flavors (four Weyl or Majorana fermions 2 ) 
transforming according to the adjoint representation of the gauge group. The global 
quantum symmetry group is £77(4). A Majorana mass term reduces this to SO (4). 



The lattice results that we report here are a continuation of our earlier work la. Il9l|. 
We will present results with spatial volumes of 8 3 , 12 3 , 16 3 and 24 3 , at several values 
of the quark mass. We are able to see trends in the finite size effects that are quite 
intriguing. 

In Section [IT] we start with a brief summary of the analytical predictions, making 
use of the old and new approaches. In addition we discuss the effect of introducing 
a non-zero fermion mass in a gauge theory that is otherwise conformal. Section II I II 
reviews our lattice simulation results. This is followed by an interpretation in Section 
IIVI Some experiments with clover fermions are summarized in Section |V] and we 
find that this improvement should allow us to explore the e and 6 regimes in a future 
work. We conclude with a discussion in Section |VT1 



II. THEORETICAL CONSIDERATIONS 



We discuss here SU(N) gauge theory with Nf Dirac fermions in the adjoint rep- 
resentation, and the critical number of flavors NJ below which scale invariance is 
broken. Because the quarks and gluons are in the same representation, it is reason- 
able to assume that Ny is independent of the number of colors N c (this is certainly 
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received recent attention 

2 Since we add a mass term in the lattice formulation, it is the Majorana description which is more 
appropriate. 
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true in perturbation theory). 



A. Truncated Schwinger-Dyson 

The first analysis of the phase diagram with fermions in higher dimensional rep- 
resentations used truncated Schwinger-Dyson (SD) equations, with fermions in the 
two- index symmetric or antisymmetric representation [6]. For N c — 2, the two-index 
symmetric representation is the adjoint representation. The generalization of the SD 
approach to any representation was carried out in [sj], yielding for the adjoint the 
conformal window 

2.075 < Nf D < — = 2.75. (2.1) 



In terms of Weyl fermions the window becomes 4.15 J$ Nyy^ D < 5.5 where the upper 
limit is the number of flavors above which asymptotic freedom is lost. The lower 
limit corresponds to the point when the SD equation can no longer be trusted and 
the anomalous dimension of the mass term is close to unity. Thus a theory of five 
Weyl fermions in the adjoint representation would appear to be in the conformal 
window, but one is uncertain what really occurs for the case of four Weyl fermions - 
equivalent to Nf = 2 Dirac fermions. Our lattice study seeks to address this question; 
however, we first describe another analytical estimate. We mention in passing that the 



approach developed in 32| provides no useful constraint for any theory with fermions 



in higher dimensional representations as shown in 33|, |34| 



B. All-orders beta function and anomalous dimension 



Specializing the recently conjectured "all-orders beta function" 11] to fermions in 
the adjoint representation, the /3-function reads: 

m = * ft-fw) , ,, 2) 



1 - &N C (l + f ) 



where j3' = N c (l — Nf) and j3o = ^(H — 4iV/) is the one-loop coefficient. 

The all-orders beta function satisfies a number of consistency checks, (i) The 
(exact) super- Yang-Mills result is recovered for Nf = 1/2. (ii) It compares well 
with the running of the Yang-Mills coupling constant as determined by lattice gauge 
theory, (iii) It provides predictions consistent with the SD approach for a critical value 
of 7 C = 1. (iv) The conformal window matches the one obtained by a conjectured 
dual gauge theory. 
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Item (iy) relates to recent exact solutions of the 't Hooft anomaly matching con 



ditions [35j, |36j. Further developments appeared in [37j|. Naturally (I2.2p reduces to 
the well-known two-loop beta function one when expanding to 0(g 5 ). We give it here 
since we will compare to it in the discussion below: 

^-A^-ik^' (2 - 3) 

with (scheme independent) adjoint representation coefficients 

iV N 2 
/3o = -£(U-4N f ), p! = -^-(34- 32JV». (2.4) 

We will also make use of the anomalous dimension 7 = — dlnm/dln/i of the renor- 
malized mass m to second order 

6iW 2N e (5W e -5N,)<f 6 
7 " (4tt) 2 + 3(4tt)* + U{9 ] {Z - b) 

The all-orders beta function predicts the anomalous dimensions of the fermion 
mass at the infrared fixed point and is in this sense "physical." In llj it was argued 
that the size of the conformal window is determined by the largest value allowed for 
the anomalous dimension, 7 C : 

11 < Nf < X" ( 2 - 6 ) 



2(2 + 7 



The fixed point value is 



11 — 4JV> 



7# = _ . (2.7) 

[ 2N f V ; 

(It is interesting to note that if we take Nf = 2, then 7* = 3/4 and the lower end of 
the conformal window in (12.61) is exactly Nf = 2.) If we use the SD inspired condition 
7 C = 1 we would have 

1.83<iV / BF <^ (2.8) 
whereas the maximal conformal window is achieved in the unitarity limit 7 C = 2: 

1.375 < Nf F < H (2.9) 

What is important to notice is that independently of which of these 7 C is chosen, 
the prediction ^, [ll| is that the adjoint theory with Nf = 2 has an IRFP and an 
associated anomalous dimension 7 = 3/4. Our aim here is to scrutinize this prediction 
of the all-orders beta function using lattice techniques. 
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FIG. 1: Beta functions for different values of the number of Dirac flavors in the adjoint 
representation of the SU{2) gauge group. The black solid curve corresponds to Nf = 2.5, 
the red to Nf = 2, the dashed one is the two-loop beta function for Nf = 2 again, while 
the magenta curve corresponds to Nf = 1.5. The green curve is the beta function for super 
Yang-Mills. 



C. Large anomalous dimensions at weak coupling 

Often, in the literature, one finds plotted a cartoon of the running of the cou- 
pling constant for either conformal or nearly conformal theories. Here we provide yet 
another cartoon of this running but this time using the "physical" form of the con- 



jectured all-orders beta function ll|. In order to be explicit, we augment this with 
a simplifying ansatz for the dependence of the anomalous dimension on the coupling 
constant, Eq. (12.5ft . The advantage is that we will be able to plot what happens when 
changing the number of flavors — but this is only a "cartoon" since we do not actu- 
ally know what 7(<? 2 ) really is in the true theory. Thus, suppose we use the two-loop 
expression of the anomalous dimension together with the all-orders beta function. 
Then the beta function for various values of Nf are shown in Fig. [TJ We have also 
plotted the two-loop beta function for SU(2) with Nf = 2. Interestingly the fixed 
point is reached before the one from the two-loop beta function. (This is consistent 



with the recent lattice results obtained in 25J.) 



What one should note from Fig. [T]is that a relatively small value of g 2 /(167r 2 ) is 
obtained while 7* = 0(1). Visually, this is because there is a long renormalization 
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group trajectory that must be traversed in going from ag « weak coupling value to 
the fixed point. The curve deepens as the number of flavors is decreased from Nf = 
2.5,2,1.5, consistent with the ordering of the fixed point values 7* = 1/5,3/4,5/3. 
The general message is that one can have large values of the anomalous dimensions 
and yet have coupling constants at the IRFP which are small. 



D. The chiral condensate 

We have measured the chiral condensate through the GMOR relation: 

{m,U? = -2m,S. (2.10) 

Here £ = (ipip) is the condensate in infinite volume. The GMOR relation just follows 
from chiral symmetry breaking with a small source m q for the "scalar current" ipip. 
In the case of an IRFP £ must also vanish as m q — > 0. Different scenarios for how it 



vanishes have been discussed in 38j. The generic expectation is that in theories where 
the anomalous dimension 7 < 1 or theories where instanton effects are important such 
as the model analyzed here, S ~ m q Afj with Ajj a high energy scale characterizing 
the onset of asymptotic freedom. In contrast QCD-like theories with chiral symmetry 
breaking possess a non-vanishing condensate as m q — > in infinite volumes. 



E. Finite volume effects 

Lattice simulations are necessarily performed on a finite four-dimensional volume, 
which we will denote L 3 x T, associating T with the extent of the temporal dimension. 
If the theory possesses an IRFP, then in the chiral limit m q —> 0, large finite-size effects 
will always be present. 

It is well-established that there are three regimes possible for lattice gauge theo- 
ries with spontaneous chiral symmetry breaking; the p-regime where m n L ^> 1 and 
m^T 3> 1, the e- regime where m n L <C 1 and m n T <C 1 and the 5-regime with 
m^L <C 1 but m n T ^> 1. In the small volume 5 or e regime the chiral condensate will 
typically scale to zero linearly with quark mass in a manner similar to that expected 
in a theory with a IRFP. Thus we must be particularly careful in interpreting our 
lattice results on small boxes in order to distinguish the two scenarios. 

Indeed, it is not clear that this categorization will prove useful in a theory that 
has an IRFP, where there is no spontaneous chiral symmetry breaking. The difficulty 
is that the expansion parameters and mode decoupling arguments rely heavily on 
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f n 7^ in the chiral, infinite volume limit. This is not true for the theory with an 
IRFP. 



III. LATTICE ANALYSIS 



For details of our lattice action and simulation algorithm we refer the reader to |l9j ]. 
Suffice it to say that we have employed unimproved Wilson fermions in the adjoint 
representation and a simple Wilson plaquette action for the gauge field and generated 
configurations using the usual HMC algorithm. We now turn to the extraction of 
accurate estimates of the meson and quark masses and the pion decay constant f n . 



A. Current quark mass extraction 

The fermion mass m q is obtained from a fit to 

n u\ 9 t GAp{t) 2Z m Z P m q _ n^+^T^ to i\ 

GpcAcit) = -p, — 777- ~ ~ " = 2m PC AC, < t < T. (3.1) 

Lrpp{t) Zj A 

Here, T is the number of sites in the temporal direction, mpcAc is the bare PCAC 
mass and m q is the renormalized current quark mass. In this work we do not determine 
the renormalization constants Z m , Zp, Z A ; however they are expected to be 0(1) and 
we will suppress them in much of what follows. The two Green's functions involved 
in (13.11) are: 



G%(t) = J d 3 x <P a (t,x)P fe (0,O)>, G%(t) = J d 3 x (^(t,x)P b (0,O)>, (3.2) 

where P a = ipj^t""^ and Aq = ^7o7st a ^, with t a £ {o" + , o" - , cr 3 }. For brevity we sup- 
press the isospin indices a, h and leave it as implied that the nonvanishing components 
G + ~ (i.e. the ones without disconnected diagrams) of the Green's functions are used 
in the measurements. At leading order in the expansion of states, 



Gpp(t) ~ d t GAp(t) ~ cosh 



T 

I — — t 



< t < T. (3.3) 



This is why a constant is expected in (13. ip . The integral over x in Eq. (I3.2p projects 
onto zero momentum states. In practice we approximate dt ~ , the symmetric 
difference operator. 
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B. Meson masses and decay constants 

Next we describe how and f n are measured. Referring to the correlation 
function Gpp(t), we work in the leading exponential approximation: 

Gppit) = Cp P cosh (m* (!-«)) ( 3 - 4 ) 

On the other hand, using the resolution of the identity in terms of states, 

G PP (t) = -j- |<0|P(0,0)|7r,<z = 0)| 2 2e-^ T / 2 cosh Cm n -tY\ (3.5) 

neglecting excited state contributions. The matrix element (0|P(0, 0) |7r, q = 0) is well 
known: 

(0|P(0,0)Kq = 0)= m ' /ff (3.6) 



Thus we obtain: 



Cj 



1/2 



U = f^j 2Z m Zpm 9 e^ T / 4 (3.7) 
This is then combined with (13.11) to obtain /^ arc = f^/Z^. 

C. Results at /3 = 2.05 

First we discuss results from simulations at (3 = 2.05 on a L 3 x 32 lattice with 
periodic boundary conditions imposed in all directions. For L — 8, 12, 16 a total of 
10500 HMC trajectories were generated. For L = 24, a total of approximately 5000 
HMC trajectories were generated. In all cases, the first 200 were discarded as ther- 
malization and observables were obtained by averaging from every tenth trajectory 
of the remaining ensemble. Errors were corrected for autocorrelations in the data. 
We performed nonlinear fits to (13.31) . We estimated the statistical uncertainty by the 
jackknife method, fitting repeatedly with a block of data removed. Jackknife block 
sizes of 10, 25, 50, 100 and 200 all gave consistent results, including error estimates. 

The fit results depend significantly on the range oft that is included, due to excited 
state contamination. The variable tfi rst determines the first timeslice that is included. 
An example of the fit variation is given in Fig. [2j It can be seen that there is a plateau 
that is reached as ifi rst is increased. In practice we fit to a constant plus exponential 
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FIG. 2: Estimates of m n a for (3 = 2.05, 8 3 x 32 lattice as the beginning of the fit range, 
tfirst, is varied. 



L 


ma 


m n a 


m p a 


fna 


m q a 


Ra 2 


8 


-1.29 


0.5686(9) 


0.6010(6) 


0.619(7) 


0.1078(5) 


11.08(10) 


8 


-1.30 


0.3680(10) 


0.3952(12) 


0.668(10) 


0.0685(3) 


13.40(15) 


8 


-1.31 


0.1433(6) 


0.1530(9) 


0.738(9) 


0.0253(2) 


19.55(8) 



TABLE I: /3 = 2.05, L 3 x 32 PBC lattice, with unimproved Wilson quarks. 

decay with tfi rst , and use the constant with fitting error as our estimate of a given 
observable. 

Results are shown in Tables [I] (mass variation) and [Til (volume variation). The 
quantity R will be discussed in a subsequent section. It can be seen that /„- decreases 
significantly with volume, whereas the behavior of the pion and rho masses are more 
complicated. Further numerical analysis is presented in Table IHIt which shows that 
a significant splitting of the rho and pion occurs for large enough volume and small 
enough quark mass. This table also shows that we are far from the heavy quark limit 
where m n « 2m q would hold. 

D. Results at j3 = 2.5 

Next we discuss results for f3 = 2.5. Here we have also allowed for a larger time 
extent, T = 64, in order to account for what might be a finer lattice spacing. The 
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L 


ma 


m q a 


m n a 


m p a 


fna 


Ra 2 


8 


-1.31 


0.0253(2) 


0.1433(6) 


0.1530(9) 


0.738(9) 


19.55(8) 


12 


-1.31 


0.015236(63) 


0.1215(17) 


0.1547(24) 


0.4598(29) 


13.83(14) 


16 


-1.31 


0.01214(16) 


0.1075(15) 


0.1531(25) 


0.406(10) 


13.854(76) 


24 


-1.31 


0.00800(11) 


0.1254(42) 


0.1770(59) 


0.1743(46) 


8.25(25) 



TABLE II: Quantities of interest for the (3 = 2.05, L 3 x 32 PBC lattice, with unimproved 
Wilson fermions. 



L 


m q a 


m n /m q 


K - m^/mn 


8 


0.0253(2) 


5.664(51) 


0.0677(77) 


12 


0.015236(63) 


7.97(12) 


0.273(27) 


16 


0.01214(16) 


8.86(17) 


0.424(31) 


24 


0.00800(11) 


15.68(57) 


0.411(67) 



TABLE III: Pion mass and rho-pion splitting enhancement as m q and 1/L are decreased, 
for the (3 = 2.05, L 3 x 32 PBC lattice, with unimproved Wilson fermions. 

fits versus tfj rst are similar to Fig. [2j The method of simulation, sampling and fits are 
the same as for (3 = 2.05. Results are given in Table HVl 

IV. INTERPRETATION 

In this section we characterize the numerical results. First consider the quantity 
R = (m 7T f n /m q ) 2 = H/m q . By our earlier arguments, for a theory with an IRFP, 
this should approach a constant in the chiral limit m q — > 0. In a QCD-like theory, 



L 


ma 


m^a 


m p a 


fna 


m q a 


Ra 2 


8 


-1.1 


0.13625(7) 


0.14531(5) 


1.039(12) 


0.03834(3) 


13.621(15) 


12 


-1.1 


0.12260(7) 


0.13537(15) 


0.593(11) 


0.02900(11) 


6.091(6) 


16 


-1.1 


0.1204(2) 


0.1251(7) 


0.405(5) 


0.0284(4) 


2.957(13) 


24 


-1.1 


0.1344(12) 


0.1497(6) 


0.242(2) 


0.0266(3) 


1.49(3) 



TABLE IV: Quantities of interest for the f3 = 2.5, L 3 x 64 PBC lattice, with unimproved 
Wilson fermions. 
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FIG. 3: Nearly conformal flow for a theory that would have an IRFP when m = 0. The 
dynamical scale "Aqcd" is generated below the mass scale m. 

the order of limits (chiral versus thermodynamic) matters. If L — > oo before taking 
the chiral limit, then R would diverge inversely with the quark mass m q . In the 8- 
or e-regime one has instead S ~ m q , which would lead to a finite result for R in 
the chiral limit. Our values of m n L at the pseudo-critical values of bare mass ma 
range from 1.1 to 3.0, whereas in the 5- or e-regime one has m n L <C 1. Thus on our 
larger lattices we are certainly outside of these small volume regimes. In Tables HT1 
and II VI we observe R decreasing as L increases. This stands in stark contrast to what 
would happen in a QCD-like theory that is outside of the small volume regimes. For 
this reason we find that our data favors the IRFP interpretation, as far as the chiral 
condensate is concerned. In fact, we find that i? ~ 1/L 2 , as a result of the observed 
scaling f n ~ 1/L. This seems to be evidence for a vanishing condensate at large 
L, consistent with the existence of an IRFP. However, we cannot rule out a small 
but nonvanishing infinite volume condensate; i.e., it is possible that we could just be 
seeing a decreasing finite volume effect that is much larger than the infinite volume 
piece for the lattices we are studying. 

The vanishing of the decay constant /„. with increasing lattice volume contrasts 
starkly with the approximate volume independence of the pion (and rho) masses 
shown in Tables [Til and IIVI The origin of this behavior is difficult to understand. 
However, we can assume that at distances scales longer than the inverse quark mass 
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the fermions effectively decouple from the dynamics leaving an IR theory which be- 
haves like quenched QCD but with a light scale 

Aqcd ~ m q e-^l£ (4.1) 



26 at 



(see Figure. [3]). These light gluonic states have been seen in simulations 
energies below that of the corresponding pion and rho states. Notice also that while 
the pion and rho masses do not scale with the lattice volume they are not simply the 
sum of the two constituent quark masses m q , as emphasized in Table IIIII Instead 
it is best to think of the pion as being composed of two "dressed" quarks where the 
dressing represents the effects of these light gluonic degrees of freedom. A useful 
analogy would the 0-system in QCD composed of two strange quarks. 

We can parameterize the dependence of the pion mass on the quark masses in a 
phenomenological way as 

m n = cm]-*L-* + C(aAuv)- ( 4 - 2 ) 

However, this formula only explains the increase in meson masses in going from L = 16 
to L = 24 provided the 0(aA^y) effects become larger in this limit, presumably due to 
larger renormalizations as the number of degrees of freedom is increased. Also, in the 
presence of an IRFP there will be large cutoff effects because one is explicitly breaking 
conformality, an essential symmetry of the theory we are trying to study. The pion 
is no longer a pseudo-Nambu-Goldstone boson, due to the absence of spontaneous 
chiral symmetry breaking, hence its mass is quite sensitive to this explicit breaking of 
scale invariance. For this reason one expects large 0(ah^y) effects, which is precisely 
what we observe. To be consistent with our numerical results the (positive) exponent 
S that appears in (14.21) should be small. 

Compatibility with the GMOR relation then implies that 

U = dm\L & - x + 0(aAl v ) (4.3) 

which has the merit of guaranteeing that f n vanishes both with the lattice volume 
and also as m q — > 0. Because f n is taken from a ratio of the lattice derivative of a 
correlation function to another correlation function, it is possible that large 0(aA^jy) 
cutoff effects may be absent, due to cancellations. In fact, this would explain our 
lattice data in the tables above, where f n is seen to decrease both with m q and 1/L; 
a large constant term C(aAu V ) does not seem to be present. Thus it may be that the 
leading lattice spacing correction to f n is actually 0(am 2 q ). 

We note that our (3 = 2.5 data is roughly consistent with Eqs. H4.2[) and (14. 3p . 
while the f3 = 2.05 is less so. For either value of /3, the increase in the meson masses 
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in going from L = 16 to L = 24 is quite strange, and seems on its face to be at odds 
with variational arguments. 3 However, topological features could play a nontrivial 
role in such finite volume considerations, in a way that might resolve the apparent 
paradox. 

V. VALENCE CLOVER ON UNIMPROVED WILSON SEA 

In order to proceed to light quark masses and move into the 5- and e-regimes 
we have experimented with simulations that utilize a clover-improved Wilson-Dirac 
propagator. We have computed the pion mass on the same unimproved dynamical 
Wilson configurations described in Section IHIl Setting the coefficient of the clover 
term csw — ®(X), we have been able to achieve m n a ^ 0.05 by tuning the valence 
quark mass. This indicates that dynamical clover fermions would allow for an explo- 
ration of the 5- and e-regimes, something that we plan to do in a future work. 

VI. DISCUSSION 

We have presented results for the low lying meson masses, decay constants and 
chiral condensate from simulations of the minimal walking technicolor theory corre- 
sponding to two flavors of adjoint Dirac fermions in SU(2) gauge theory. Data at 
two couplings j3 = 2.05 and (3 = 2.5 and a range of lattice volumes L 3 x 32(64), 
with L = 8, 12, 16, 24 were shown. Unimproved Wilson fermions and Wilson glue are 
used and ensembles of 0(10000) configurations accumulated at each set of parameter 
values. 

We have shown how the GMOR relation may be used to compute the chiral conden- 
sate and discussed the relationship between the condensate measured on the lattice 
and its continuum cousin. Our results are consistent with the vanishing of the con- 
densate in the infinite volume limit and hence the existence of an infrared fixed point. 
The dependence of f n and the pion and rho masses on both the quark mass and 
lattice volume are shown to also support the presence of such a fixed point, though 
other interpretations are possible. 



3 We thank R. Brower for raising this point. The essence of the argument is that on doubling the 
lattice size, the original pion wavefunction can be periodically extended. But one would expect 
that it is no longer the minimum energy eigenstate in the pseudo-scalar channel, since new basis 
states for a variational analysis are allowed on the larger lattice. It follows that the pion mass in 
the larger volume will be lower. 
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